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ABSTRACT 



By cross-correlating large samples of galaxy clusters with publicly available radio source catalogs, we construct the volume-averaged 
radio luminosity function (RLF) in clusters of galaxies, and investigate its dependence on cluster redshift and mass. In addition, we 
determine the correlation between the cluster mass and the radio luminosity of the brightest source within 50 kpc from the cluster 
center. We use two cluster samples: the optically selected maxBCG cluster catalog and a composite sample of X-ray selected clusters. 
The radio data come from the VLA NVSS and FIRST surveys. We use scaling relations to estimate cluster masses and radii to get 
robust estimates of cluster volumes. We determine the projected radial distribution of sources, for which we find no dependence on 
luminosity or cluster mass. Background and foreground sources are statistically accounted for, and we account for confusion of radio 
sources by adaptively degrading the resolution of the radio source surveys. We determine the redshift evolution of the RLF under the 
assumption that its overall shape does not change with redshift. Our results are consistent with a pure luminosity evolution of the 
RLF in the range 0.1 < z < 0.3 from the optical cluster sample. The X-ray sample extends to higher redshift and yields results also 
consistent with a pure luminosity evolution. We find no direct evidence of a dependence of the RLF on cluster mass from the present 
data, although the data are consistent with the most luminous sources only being found in high-mass systems. 

Key words. Galaxies: clusters: general - Radio continuum: galaxies - Galaxies: active - Galaxies: evolution 



1. Introduction 

The study of radio sources inside galaxy clusters is becom- 
ing increasingly import ant in the context of large-area Sunyaev- 
Zel'd ovich effect (SZE. ISunvaev & Zeldovich 1972; Birkinshaw 



19991) cluster surveys (e.g. ISchwan et all 120031: iHo et al.l 120091: 
Vand erlinde et alJ[2010t lMuchovei et al.ll2010t) . especially in the 



case of single-frequency SZE surveys and where the resolution 
is matched with the typical cluster size. Powerful radio sources 
are generally associated with early type galaxies, and as the latter 
preferentially reside in clusters, it is expected that radio sources 
can alter the SZE decrement in both position an depth. 

Apart from contamination of SZE surveys, a study of the 
radio-loud population of AGN inside galaxy clusters sheds light 
on the interaction between c luster cooling flows and AGN heat- 
ing/feedback scenarios (e.g. Ma rtini et al.ll2009l) . Analyzing the 
radio properties of galaxy clusters can provide direct evidence 
of heating of the intra-cluster medium (ICM) by AGN through 
various stages of cluster formation (e.g. lHart et al.ll2009l) . 

Three pieces of information are crucial to assess the radio 
source contamination of SZE observations: (1) the distribution 
of sources as a function of projected distance from the cluster 
center; (2) the brightness distribution of sources and its redshift 
evolution; and (3) the spectral energy distributions (SEDs) of the 
sources. A convenient way of dealing with the first two points is 
to construct a radio luminosity function (RLF), i.e. the number 
density of sources as a function of luminosity, averaged over the 
cluster volume. This is the objective of the present paper. 
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iLedlow & Owen] {l996) first constructed the RLF for radio 
galaxies in clusters using a 20 cm VLA survey of Abe l l clus- 
ters. Similar studies were carried out by Pedd v & Yunl (12 004) 
for sev en nearby galaxy clusters, and by Massardi & De Zotti 
(12004]) fo r a much larger sam ple of 951 Abell clusters at low 
redshifts. Sto cke et all (1 19991) used a sample of 19 X-ray se- 
lected clusters at 0.2 < z < 0.8 to constrain the evolution 
of radio galaxie s , and comparing their results with those of 
ILedlow & Owen! ([1996), they found no evidence for an evolu- 
tion of the radio-loud AGN po pulation. Similar conclusions were 
drawn by ICoble et all (120071) using a sample of massive clus- 
ters detected through t he SZE decrement at 30 GHz. Although 
Branchesi et al.l d2006l) found indications of redshift evolution 
in both slope and amplitude of the cluster RLF when comparing 
VLA observations of a sample of 18 X-ray selected clusters with 
the results of Stocke et al., no definitive conclusions were drawn 
on the RLF redshift evolution. 



Recent high-resolution X-ray imaging of clusters have pro- 
vided evidence that t he AGN fraction inside clusters rapidl y 
evolves with redshift (lEastman et alJl2007t iMartini et al.l l2009). 
These results are affected by small number statistics, since the 
AGN fraction in any single cluster is very small and can also 
vary with other cluster p roperties, such as the velocity disper- 
sion (Sivakoff et al. 2008). The volume-averaged radio luminos- 
ity function is therefore expected to be a more robust indicator 
of redshift evolution. Correla ting X-ray selected cl usters with 
optical and IR selected AGN, iGalametz et alJ ((2009) confirmed 
that the AGN excess near cluster centers increases with redshift. 
A parametric modeling of this positive evolution, which can be 
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used to predict the radio AGN contamination in SZE cluster sur- 
ve ys, was not considere d in these studies. 

iLin & Mohri (120071 hereafter LM07) used a sample of 573 
X-ray selected clusters to estimate masses and virial radii, yield- 
ing a more physically meaningful cluster volume for the RLF 
than obtained when using a constant cluster radius. We adopt 
that method in this paper. LM07 found a narrow radial distribu- 
tion of radio sources, which we reproduce in this work. 

This study focuses on the redshift evolution of the RLF and 
its decoupling from a possible dependence on cluster mass. We 
use publicly available radio source catalogs at 1 .4 GHz to iden- 
tify radio sources associated with clust ers of galaxies in a large 
sample of optically selected clusters dKoeste r et al. 2007) and a 
composite sample of ROSAT selected X-ray clusters. Studies of 
radio point sourc es in clusters at higher frequencies were car- 
ried out by, e.g ., ICoorav et ail (1 19981) . ICoble et al.1 (120071) and 
ILin et all (120091) . and the extrapolation of the RLF from 1 .4 GHz 
to frequencies relevant for SZ E surveys was investigated by, e.g., 
iMassardi & De Zottil (120041) and by LM07. These two aspects 
are not pursued in this work. 

This paper is organized as follows: in ^2] we describe the 
samples of radio sources and cluster of galaxies. We discuss the 
estimation of cluster masses from observable properties for clus- 
ters selected from optical and X-ray cluster catalogs, and de- 
scribe how radio luminosities are derived from flux densities. 
The radial density profile of radio sources associated with clus- 
ters of galaxies is discussed in $3] In $4] we investigate the cor- 
relation between cluster mass and the luminosity of the brightest 
radio source in the central region of the cluster. The main part 
of the paper is where the 1 .4 GHz volume-averaged radio 
luminosity in galaxy clusters is derived from our sample. The 
method of computing the RLF is described, and the statistical 
treatment of sources not associated with the clusters is discussed. 
Confusion of sources is discussed, and the redshift evolution of 
the RLF is derived. We summarize our main results and offer our 
conclusions in $6] 

Where not otherwise noted, we use the cosmological param- 
eters from the WMAP 5 -year cosmology dKomatsu et al.l [2009) 
with h = 0.705, Q. m h 2 =0.136 and Q A = 0.726. 



2. Cluster and radio source samples 

In this section we describe the galaxy cluster samples and 1.4 
GHz radio source catalogs used in this paper. 



2.1. Cluster samples 

To demonstrate a redshift evolution of the RLF in the best pos- 
sible way, we select clusters with a wide range of redshifts from 
publicly available optical and X-ray surveys. To make robust es- 
timates of the volume averaged RLF, we derive cluster masses 
and radii using scaling relations. The mass within a region whose 
mean density is 200 times the critical density of the universe at 
the cluster redshift, M?oo , has been f ound to be a good estimator 
of the virial mass re.g. fmit3l200l . and is used here together 
with the corresponding radius, r2oo, which we use to define the 
typical scale of the cluster. 



2.1 .1 . Description of the samples 

Because of its large sampl e size, we use as our main sample the 
optical maxBCG catalog dKoester et al J 1200*71) . which contains 
a total of 13,823 clusters extracted from the Sloan Digital Sky 
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Fig. 1. Distribution of X-ray selected clusters of galaxies in the 
Lx- z plane. 

Survey (SDSS). This is by far the largest homogeneous publicly 
available cluster sample, and its redshift range 0.1 < z < 0.3, 
corresponding to approximately 2.5 xlO 9 years of cosmic time, is 
large enough to allow an investigation of the redshift dependence 
of the RLF. 

To allow a comparison with X-ray selected clusters, we 
have gathered in one large sample most of the published clus- 
ter detections from the ROSAT mission. Data products from this 
satellite remain the reference to date since only few new mas- 
sive systems have been reported from more recent observato- 
ries. Our source list includes the ROSAT all-sky surve y cata - 
logs from NORAS and REFLEX dBohringer et al.ll2000l BOOl . 
and some serendipitous catalogs ext racted from the archiv es of 
pointed o bservations: the 160 deg 2 dVikhlinin et al.lll998l). the 
400deg 2 dBurenin et all 120071) and WARPSI/II dBurenin et all 
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Horner et al.l 120081) catalogs. We also added the com- 



plete sample of 12 high redshift very luminous systems of 
lEbeling et alj d2007l) . In total this yields 1177 X-ray selected 
clusters. The number of clusters from the respective catalogs and 
the corresponding flux limits are given in Table [T] Mixing these 
catalogs enables us to sample a large fraction of the Lx —z plane, 
as can be seen in Fig.Q] 

2.1.2. Cluster masses and radii 

For the maxBCG sample, we esti mate halo masses from 
the scaled richness parameter n^ 00 . iKoester et alJ d2007l) and 
Beck er et al .1 (120071) find that this measure correlates well with 
the gal axy velocity dispersion, and thus also with halo mass. 
Indeed, Andre on & Hurnl d2010h found that optical richness and 
X-ray luminosity perform similarly well in predicting cluster 
masses. The relation between nf. 2 , 00 and M200 has been exten- 



sively studied, both using weak l ensing (Johnston et al. 2007; 
iMandelbaum et alj2 008: She ldon et alj2009HRozo et alj20 09a) 
and comparing to X-ray luminosities (Rvkof f et alj l2008al t : 
iRozo et al . 20093). 

In this paper, we use the X-ray luminosity-op tical richness 
relation (L x - n s20 °) found bv lRvkoff et ail d2008bl) to derive the 
expected X-ray luminosities as 



LxiRwo) = e" 
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Table 1. Our composite X-ray cluster sample. The third column is the effective flux limit, in the [0.5-2] keV energy band, estimated 
from our recomputed source parameters. 



Catalog 


Number of clusters 


Limiting flux (erg s 1 cm 2 ) 


REFLEX 


447 


2.1e-12 


NORAS 


371 


1.2e-12 


160deg2 


221 


8.0e-14 


400deg2 


242 


1.4e-13 


WARPS 1+2 


124 


8.0e-14 


MACS(z>0.5) 


12 


7.0e-13 



where a = 3.90 + 0.04, yS = 1.85 ± 0.05 and y = 6.0 ± 0.8. We 
scale the derived X-ray luminosities to mass using a self-similar 
redshift evolution, in the same way as described below for the 
X-ray sample. 

For the X-ray sample, the heterogeneity of the data makes 
it difficult to straightforwardly estimate the cluster masses in 
a uniform way. The most basic observable, available for all 
the sub-samples although in different bands {[b m i n ,b max \), is 
th e total cluster X-ray flux. We thus apply the scaling relation 
of lArnaud & Evrardl (11999b . linking the average gas tempera- 
ture T to the bolometric X-ray luminosity L\, Q \, with a self- 
similar redshift evolution to define a unique mapping between 
(z,Fx,b m i n ,b max ) and (Li,„i,T). For th is, we make use of the APEC 
spectral model (Smith et al. 2001) with heavy element abun- 
dances set to 0.3 times the solar values. The derived cluster tem- 
peratures are translated into M200 following the scaling law of 
Arnaud et al. (2005), again assuming self-similar evolution. 

From simple geometric considerations and the definition of 
over-density with respect to the critical density, it follows that 
r2oo is related to M200 by 



M 200 = (4^/3)4,0 200p c (z) 



(2) 



where p c (z) is the critical density of the universe at redshift z. 
r2oo is typically in the range 1-3 Mpc for the clusters in our 
samples. 

2.1.3. Cluster sample selection criteria 

In order to limit our survey to massive systems, we exclude all 
clusters less massive than 5 x 10 13 M o from our sample. We also 
exclude low-redshift systems with z < 0.1, in order to have the 
same redshift cutoff in both samples. In addition, we define a 
low-redshift (0.05 < z < 0.12) sub-sample of the X-ray sam- 
ple, with the same mass cut (see Tabled, to investigate the mass 
dependence of the RLF ( 95.3.31 ). The lower redshift limit of the 
sub-sample is chosen so as to avoid excessive overlap of low- 
redshift cluster fields. The upper redshift limit is chosen as high 
as possible without including redshifts where the mass limit of 
the flux-limited X-ray cluster sample is greater than our mass 
limit of 5 x 10 13 M G , in order to avoid a bias in the the determi- 
nation of a possible mass dependence in the RLF. The fact that 
the low-redshift sub-sample has a slight overlap with the high- 
redshift X-ray sample is of little consequence as the two samples 
are never compared directly with one another. 

Not all cluster fields in our samples are covered by the NVSS 
and FIRST radio surveys that are used for this study. Apart from 
taking this into account, we note that due to dynamic range limi- 
tations of the radio interferometric data, in the FIRST and NVSS 
catalogs some regions on the sky around strong radio sources 
are plagued by abnormally high or low source counts, often with 



catalog entries not corresponding to real objects{3 We identify 
these regions by counting the number of sources in a circular 
region with radius 1 around each cluster center, and excluding 
fields where the source counts exceed or fall below the average 
counts in cluster fields by more than three standard deviations^ 
Note that using the average counts of FIRST and NVSS would 
have underestimates the expected counts in cluster fields, where 
a local surface over-density is expected. Counting sources within 
one degree of maxBCG clusters (cf. @, we find the average 
density of FIRST sources in such a region to be 9% higher than 
the catalog average of 90 sources per square degree. In the case 
of NVSS, only about 4% of the cluster fields are affected by this 
cut (determined from the difference of rows 4 and 5 of Table [2] 
in the maxBCG case). 

We also exclude pairs (or multiplets) of clusters where the 
sum of the radii (r2oo projected on the sky plane) is greater 
than the angular separation between the cluster centers (row 2 
ofTable0. 

Table|2]summarizes the selection criteria applied to the clus- 
ter sample, and lists how many clusters in each of the final sam- 
ples are covered by the FIRST and NVSS catalogs. 

2.2. Radio source samples 
2.2.1 . Description of the samples 

To search f or radio sources associated with clus ters, we use 
the FIRST (iBecker et alJ [19951 IWhite etail fl99l and NVSS 
dCondon et alJ 1 19981) 1.4 GHz radio continuum surveys. The 
most relevant properties of the NVSS and FIRST surveys for 
the present work are summarized in Table [3] We make direct use 
of publicly available source catalogs from each survey. While 
the NVSS catalog has the advantage of greater sky coverage (cf. 
Table [3]), thereby covering more of our cluster fields, it has the 
disadvantage of a poorer resolution, resulting in increased source 
confusion. 

We search for FIRST and NVSS sources within a projected 
radius corresponding to 3r2oo around each cluster center, with 
r2oo computed from Eq. ©. In order to facilitate a robust com- 
parison of NVSS and FIRST sources, we exclude all sources 
fainter than 5 mJy, which is well above the completeness limit 
of both surveys. 

In principle the use of the FIRST catalog alone would be 
enough for the purpose of this paper. We make use of the NVSS 
catalog at all stages of the analysis because it serves as a test 
that our results obtained with FIRST are robust. In particular, the 



1 Typical examples are the fields centered on the quasars 3C273 and 
3C295 

2 The 1° radius corresponds to approximately ten times r 2 oo at the 
median redshift and mass of the maxBCG sample, and is chosen so as 
to include all cluster sources and to have a sufficient number of back- 
ground and foreground sources to allow for robust statistics 



4 



M. W. Sommer et al.: Galaxy cluster radio luminosity function 



Table 2. Cluster sample selection criteria. The selection is cumulative in the sense that in rows 3-4, each row assumes all the 
conditions of the previous rows. Rows 5 and 6 assume all the conditions in rows 2-4. The low-z X-ray sample is used only to 
constrain the mass dependence of the RLF ( 35.3.31 l. Where not explicitly stated otherwise, the high-redshift samples are used in this 
paper. 





Main sample 


maxBCG 


X- 


ray 


1 


clusters in main sample 


13823 


1177 


2 


clusters with sufficient separation 


12846 


1121 




Sub-sample 




high-z 


low-z 




Redshift range 


0.1 <z<0.3 


0.1 <z< 1.26 


0.05 < z < 0.12 


3 


Clusters within redshift range 


12846 


690 


292 


4 


Clusters with M > 5 x 10 13 M o 


12522 


674 


275 


5 


clusters with NVSS coverage^ 


12475 


596 


218 


6 


clusters with FIRST coverage" 


11812 


273 


75 



" Excluding regions with abnormally high or abnormally low source counts (see text) 



Table 3. Properties of the FIRST and NVSS radio continuum 
surveys 





FIRST 


NVSS 


effective resolution 


5" 


45" 


completeness limit 


1 mJy^ 


~2.5 mly 


positional uncertainty^ 


< 0.5' 7 


< 1" 


positional uncertainty^ 


1" 


~7" 


total number of sources 


8 16,33 Vl 


1,810,672^ 


sources per square degree 


-90 ~ 


-53 


area covered (deg 2 ) 


9055 


33885 



" catalog detection threshold 

* for the brightest sources 

c at the completeness limit/detection threshold 

d as of July 16, 2008 

' as of February, 2004 

NVSS data are important for testing our method of adaptively 
accounting for confusion, as described in 35.21 Because the two 
surveys have significant overlap, we do not need to "mix" them 
in the sense of combining source counts from the two catalogs 
for any particular sample. 



avoid counting cluster fields where the bin luminosity corre- 
sponds to a flux below our chosen threshold of 5 mJy. z cut is 
determined from Eq. d3) and only affects the lowest luminosities 
considered in the analysis. 

3. Radial density distribution of radio sources 

As we have no redshift information for individual radio contin- 
uum sources, the over-density of sources toward clusters must be 
quantified statistically. We construct the stacked radial profile of 
radio sources around cluster centers. Following LM07, we take 
r/r2oo as the radial coordinate. The resulting radial profile in- 
cludes both cluster galaxies and the field (background and fore- 
ground) population. Although it cannot be decided whether indi- 
vidual sources are cluster members, cluster and field sources can 
be separated statistically, since the field population has a radially 
constant contribution to the radial profile. Besides determining 
the stacked radial density distribution of all radio sources, we 
also investigate whether the distribution depends on source lu- 
minosity and host galaxy cluster mass. 

The radial profiles are used when constructing the luminos- 
ity function by de-projection into volume number density as de- 
scribed in 35.1.31 



2.2.2. Radio source luminosities 

To compute radio luminosities, we assume isotropic emission, 
which is not correct for any individual galaxy, but is well jus- 
tified when averaging over the entire sample. The luminosity 
(power) of a radio source is 

il.4GHz = (47rP?) S 1.4GHz "~~T » (3) 

(1+z) 

where S i.4ghz is the angular integrated flux density taken from 
the VLA catalog (FIRST or NVSS), D L is the cosmological lu- 
minosity distance and < K{z) is the ^-correction. 

The radio sources are modeled with continuum spectra of the 
form S v oc v~ a , where a is the spectral slope. When computing 
the luminosity from flux, we have to account for the fact that due 
to the redshift, the observed flux corresponds to a rest frequency 
higher than 1.4 GHz. For the computation of this so-called k- 
correction, < K(z), we assi gn a spectral slope of a — 0.72 to all 
sources, as determined bv lCoble et al.l (2007) in the range 1 .4 — 
30 GHz. The effect of this correction is about 10% for sources at 
z = 0.2 and 30% at z = 0.6. 

When binning the data by luminosity to compute the RLF, 
we determine a redshift limit, z cut , for each luminosity bin to 



3.1. Radial model 

The angular offset of radio sources with respect to the cluster 
center (defined as the position of the brightest cluster galaxy 
(BCG) in the maxBCG sample) is translated into projected phys- 
ical distance at the cluster redshift. The physical radial distances 
are stacked for all cluster fields and binned by radius. 

We fit the radial profile of radio sources in each cluster sam- 
ple by a parametric model. In general, all our radial profiles are 
well fit by a model of the form 

/ f 2 r^ i ( 2 

^) = ^h + y + -iac^-' + tf, (4) 

where is the projected surface density and ^ = r/r2oo- 
The fir st term of dp has the form o f an isothermal (3- 
model (Cavaliere & Fusco- Femiandfl978l) . the second term is 
a Gaussian, corrected for the area in an annulus, and the third 
term, </<y, is the constant field density. 

The isothermal /J-model is described by three parameters: the 
(peak) normalization \frp, the scale radius and the power law 
index /?. The /J-model is chosen because it provides a good fit to 
our radially binned data outside r/r2oo — 0.006. 
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In our fractional radial coordinates, the constant field density, 
t//f, is a function of the redshift distribution of the galaxy cluster 
sample. It is not physically meaningful as we have (incorrectly) 
placed all radio sources at the cluster redshifts, but as we are not 
interested in the field population this component may safely be 
ignored. 

The second term in Eq. © accounts for an additional 
peaked feature in the radial distribution of FIRST sources in 
the maxBCG galaxy cluster sample (Fig. [2). It is well fit by a 
Gaussian with normalization ipc and variance a 2 . The origin of 
this feature is the fact that the cluster center is taken as the po- 
sition of the BCG. As discussed further in 33.31 the underlying 
distribution derives from the extended nature of the large major- 
ity of FIRST radio sources associated with the BCG. Note that 
the central Gaussian component of if/(t;) is considered only for 
the radial distribution of FIRST sources in the maxBCG sample; 
in all other cases we set i//c = 0. 

For the purpose of constructing the RLF, the exact form of 
the fitting function is not important; although the density pro- 
file is used in constructing the RLF (as discussed in 95.1.31 ). the 
RLF is only very weakly sensitive to the parameterization of the 
source density fitting function. 

3.2. Dependence on luminosity and cluster mass 

Because the radial distribution of sources plays an important role 
in determining cluster member radio source counts within r2oo, 
it is necessary to investigate whether bright radio sources have a 
radial distribution different from the distribution of faint sources. 
We compute the luminosities of FIRST sources in maxBCG 
clusters according to Eq. (|3} by placing all radio sources at the 
cluster redshifts and binning the sample in luminosity. We deter- 
mine the radial source distribution in each luminosity bin using 
the fitting function given by Eq. Similarly, we investigate a 
possible dependence on cluster mass by binning the cluster sam- 
ple in mass. 

Because the sample has been divided and uncertainties on 
individual profiles are greater, we fix the power law parameter 
/J in Eq. © to the best-fit value f3 = 0.987 from the total fit 
(Table|4]l in order to get reliable estimates on the scale radius f . 

3.3. Results 

3.3.1. Ensemble properties 

Figure [2] shows the source density profiles in maxBCG clusters, 
fitted with the parametric form (2) and divided into its compo- 
nents. Table [4] indicates the most important parameters in the 
radial fits of radio sources associated with clusters. 

Note that the parameters {^iq,ct) of the Gaussian were fit 
only for the FIRST sources associated with maxBCG clusters, 
as this is the only combination of cluster and source catalog that 
yields an additional centrally peaked feature. The width of the 
peaked distribution (~ 3 x 1CT 3 r^oo) corresponds to 0.9" at the 
median distance of the maxBCG clusters. Although this is com- 
parable to the maximum positional uncertainty of FIRST point 
sources (~ 1") at the completeness limit, the FIRST positions 
are typically constrained to a much better accuracy, in particular 
at higher flux levels. Thus, positional uncertainties in the FIRST 
survey alone cannot explain the width of the peaked distribution. 

A visual inspection of a large fraction of ~1800 FIRST 
sources with positional matches out to 15" of the maxBCG po- 
sitions revealed that the very close matches (< 0.3") are domi- 
nated by compact radio sources, while the matches in the range 




0.001 0.010 0.100 1.000 

t = r / r 200 



Fig. 2. Radial distribution of 1 .4 GHz radio sources brighter than 
5 mJy in the maxBCG cluster sample. The normalization on 
the y-axis is arbitrary. Radio sources selected from the FIRST 
catalog (triangles with error bars) clearly indicate three compo- 
nents: a narrow central peak, fitted with a Gaussian (dotted line); 
a broader distribution of sources, fitted with a y6-profile (short- 
dashed line) and a background/foreground component (dash- 
dotted line). The sum of the components is indicated by the solid 
line. The radial profile derived from the NVSS catalog is indi- 
cated by the error bars without symbols. The fit to these data 
(red long-dashed line) does not include a Gaussian component. 



3"— 10" can be classified into three broad categories: (a) dou- 
ble sources with more or less symmetrical lobes straddling the 
BCG position, (b) triple sources with weak cores (centered on 
the BCG) and strong lobes, for which our threshold of 5 mJy 
picked up only the lobes, and (c) complex sources, usually of 
FRII or wide-angle tail type, some of which have truly spectac- 
ular radio morphologies and sizes up to 1 ' or more. From this 
inspection we estimate that out to 10" separation (~ 0.3r2oo at 
the median redshift) at least 98% of our matches correspond to 
the BCGs (with at most 10% of the sources being duplicates, 
i.e. different components or lobes of the same complex source). 
While for separations between 10 and 15" (i.e. up to 54 kpc at 
the median redshift of the clusters of z=0.227) the fraction of 
duplicates exceeds 30%, our matches can be safely related to the 
BCGs in at least 90% of the cases. We conclude that the peaked 
component in the radial distribution originates from mismatches 
in SDSS and FIRST positions caused by extended radio emis- 
sion. 

Because the BCG component can be separated from the gen- 
eral distribution of FIRST sources in the maxBCG sample, the 
dimensionless core radius, for the remaining sources is much 
larger than for the NVSS sources, where the peaked central dis- 
tribution is "hidden" in the y6-model component. The distribution 
of sources in the X-ray cluster sample is less peaked (larger £) 
because the cluster center definition is based on the gas distribu- 
tion and not the position of the BCG. 

We note that our results pertaining to FIRST sources in 
maxBCG c luster s verify the radial distribution derived by 
ICroftetalJ(l2007h for the same sample. As indicated in Fig. [2] 
our narrow binning reveals the details of the centrally peaked 
component in detail and confirm its origin as the BCG popula- 
tion. The relatively high redshifts in the optical sample (z > 0.1) 
also imply that the population of star forming galaxies in clusters 
- with a broader distribution than radio-loud AGN - has been ex- 
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Table 4. Results of j6-model fits to the radial distribution of sources associated with clusters of galaxies. The data are fitted using 
Eq. ©. 



Cluster sample 


Source catalog 


P 


e 


if/ G cr 


maxBCG 


FIRST 


0.987±0.018 


0.0115 +0.0012 


3.23±0.18 0.00274±0.00013 


maxBCG 


NVSS 


1.018+0.010 


0.00257±0.00019 




X-ray 


FIRST 


0.899±0.136 


0.0191 ±0.0099 




X-ray 


NVSS 


0.911±0.083 


0.0198 ±0.0063 





Table 5. Results of /J-model fits to the radial distribution of 
FIRST sources associated with maxBCG clusters of galaxies, 
with the radio sources binned by luminosity L. Only the param- 
eters f and cr, pertaining to the width of the profiles, are listed. 
P has been fixed to the best-fit value 0.987 from the total fit. The 
data are fitted using Eq. (0). 



Luminosity range ( cr 

log(L [W/Hz]) < 23.8 0.0116 ±0.0018 0.00249 ±0.00067 

23.8 < log(L [W/Hz]) < 24.2 0.0138 ±0.0023 0.00319 ±0.00054 

24.2 < log(L [W/Hz]) < 24.6 0.0138 ±0.0024 0.00343 ±0.00073 

log(L[W/Hz])>24.6 0.0097 ±0.0018 0.00303 ±0.00088 



Table 6. As table[5] but with the sample divided into mass bins. 



4.1. Background 

Brightest cluster galaxies (BCGs) tend to lie at the center of the 
mass distribution and are more luminous than average ellipticals. 
Moreover, their properties correlate strongly with their host clus- 
ters, a s seen from numerical simulations (e.g De Lucia & Blaizot 
2007) and optical/ near-IR observations (e.g. iLin & M ohr 2004; 
Brough etal. 2008). It has been shown that BCGs are about 
an order of magnitu de more likely than other ellipticals to 
host r adio-loud AGN dEdgdl 199 It ILin & Mohrll2004t iBest et al.l 
2007). Based on this, it is likely that the most radio-luminous 
galaxies in clusters are associated with BCGs, and in addition to 
the correlation between BCG luminosity and cluster mass one 
can expect that 1 .4 GHz radio luminosities of the central radio- 
loud AGN are also correlated with cluster mass. 



Mass range 


t 


cr 


log(M/M G ) < 13.90 
13.90 < log(M/M Q ) < 14.15 
14.15 < log(M/Af Q ) < 14.40 

log(M/M Q ) > 14.40 


0.0132 ±0.0016 
0.0139 ±0.0017 
0.0125 ±0.0011 
0.0109 ±0.0011 


0.00334 ±0.00060 
0.00344 ±0.00075 
0.00260 ±0.00053 
0.00209 ±0.00049 



eluded by our flux limit of 5 mJy. As a consequence our y6-model 
fit to the radial profile is much narrower compared to those ob- 
tained for the radio source distribution in loc al clusters (e.g. 
iReddv & Yunl2004tlMassardi & De Zottill2004 . Similar narrow 
radial profiles consisting mostly of radio-loud AGN were also 
found by LM07. 

3.3.2. Luminosity and cluster mass dependence 

From their X-ray sample of galaxy clusters, LM07 found that 
more luminous radio sources are more centrally concentrated in 
galaxy clusters. As indicated in Table [5] we are unable to re- 
produce this result with our optical maxBCG sample; although 
the highest luminosity bin has a slightly smaller scale radius f 
of the outer profile, the result is not significant. Similarly, there 
is no statistically significant indication of a mass dependence of 
the radial density from the maxBCG sample of galaxy clusters 
(Table [6j. We verify that the same is the case when using the 
NVSS sample. 



4. Mass-luminosity correlation 

In this section we investigate the correlation between cluster 
mass and the radio luminosity of the brightest radio source in the 
central region of the cluster. We compare the results with similar 
studies done with the optical luminosities of the brightest cluster 
galaxies, and discuss the implications for our goal of compar- 
ing the volume averaged radio luminosity function in clusters at 
different redshifts. 



4.2. Method 

We investigate the correlation of BCG radio luminosity with 
cluster mass by taking the luminosity of the single brightest 
source within a radius of 50 kpc from the cluster center. We reit- 
erate that the latter is actually the BCG position in the maxBCG 
sample. In the X-ray sample, the search radius allo ws for offsets 
betwe en the X-ray center and the BCG position; iKoester et al.l 
d2007l) found typical offsets of 50 /T 1 kpc by cross-correlating 
the maxBCG catalog with NORAS and REFLEX data. 

Note that this selection method leads to a different minimum 
flux level for different redshift sub-samples, since we are dis- 
regarding clusters that do not contain any radio source above a 
certain flux density. This leads to a higher mean luminosity of 
the central radio source for clusters at higher redshift, compared 
to clusters at lower redshifts within the same mass bin (derived 
from the X-ray or optical mass observables). The actual redshift 
evolution of the radio luminosities of the central BCGs is hid- 
den within this luminosity bias. Considering the large scatter in 
the luminosity of the central radio source, we do not attempt to 
model the redshift evolution from this method. 

4.3. Results 

Figure [3] shows the correlation of luminosities of FIRST sources 
with cluster mass for both our cluster samples. Modeling the 
correlation with a power-law of the form L ~ M", the best- 
fit correlation from the maxBCG sample in the redshift range 
0.1 < z < 0.2 is a = 0.31 +0.12. Although the X-ray sub-sample 
in the same redshift range is essentially too small to constrain 
the power law, the data are consistent with the maxBCG result, 
with a = 0.34 + 0.41. The low-redshift X-ray sub-sample yields 
a = 0.44 + 0.27, again consistent with the maxBCG relation. 

The found correlation between radio luminosity and cluster 
mass is consist ent with results fro m optical/near-IR observations 
of the BCGs. lLin & Mohrl (12004 compared the K-band near-IR 
luminosities for a s ample of X-ray select ed clusters and found 
Lbcg <* M° d 26±0M . lHaarsma etail (120101) reported a shallower 
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Fig. 3. Correlation between radio luminosity of the brightest 
FIRST source within a projected physical radius of 50 kpc from 
the cluster center (BCG position in the maxBCG sample) and 
cluster mass for three sub-samples of galaxy clusters. Diamonds 
indicate maxBCG clusters in the redshift range 0.1 < z < 0.2, 
and the best-fit correlation is indicated by the dotted line. This 
is consistent with the X-ray sample in the same redshift range 
(triangles). Fixing the slope, the correlation is consistent with a 
local sub-sample of X-ray clusters (0.01 < z < 0.1, squares), 
where the luminosity level is lower due to a different luminosity 
cutoff in this sub-sample. 



correlation from an X-ray selected sample; Lbcg 00 M°; °' 07 . 
Thus, luminosities in other bands show correlations with mass 
that are similar to our results. 

For the optical cluster sample, the correlation of radio lu- 
minosity and cluster mass strongly points towards increased ra- 
dio emission from the central AGN in the BCGs of more mas- 
sive clusters. ICroft et al.l(l2007l) also computed the rest-frame 1 .4 
GHz luminosities of the BCGs by cross-correlating the maxBCG 
catalog with the FIRST survey, and noted the trend of increas- 
ing radio luminosity with optical luminosity of the BCGs in the 
r-band. We do not expand on this result further, as the radio lu- 
minosity of the central brightest radio source is a poor indicator 
of the total radio luminosity in clusters. A cluster can have multi- 
ple BCGs as a result of its merging history .0 and the radial source 
density profile in Fig. [2] suggests that the radio source popula- 
tion extends beyond the cluster core region. For these reasons, 
we now turn to the construction and modeling of the volume av- 
eraged radio luminosity function in galaxy clusters. 



5.1. Method 

The luminosity function (p(L) is defined as the average number 
of radio sources per unit physical volume of a cluster and per 
logarithmic luminosity bin ("magnitude") . We use the estimated 
values of r2oo from 92.1l to define cluster volumes, expressed in 
Mpc 3 . 



5.1.1. Source counts 

Figure |4] illustrates schematically how source counts are ob- 



Observer 




200 



Fig. 4. Illustration of the method used to obtain source counts 
and volumes for the computation of the luminosity function of 
radio sources inside clusters. For each cluster the sampled vol- 
ume is the "line-of-sight cylinder". The projected number den- 
sity in this volume is converted to volume density of sources in- 
side the cluster radius (r2oo) by applying a model for the volume 
density of sources. 



tained. For each cluster, a region within a projected physical ra- 
dius r = J] r2oo, from the center is searched for radio sources 
in a given luminosity bin, using the cluster redshift to convert 
from flux to luminosity. As explained below, the search radius 
does not have a systematic effect on the RLF; however, it can 
be chosen such as to minimize the uncertainties of the latter. 
The number of sources found within the chosen radius is A^ ount . 
We need to convert this number into a de-projected number of 
sources within r2oo- For this, we need to know (i) the number of 
foreground/background sources, and (ii) the radial source den- 
sity profile. In the following we discuss these issues in turn. 



5. The radio luminosity function 

In this section, the radio luminosity function (RLF) is derived 
for each of our samples. We start by discussing how accurate 
source counts in the cluster volume are obtained Q5.U . In jj5.2l 
we describe how source confusion is accounted for. We describe 
a parametric model for the RLF in 95. 31 and also discuss how a 
redshift evolution is modeled. We give the results in 95. 41 



3 After a cluster merger, it will take considerable time before the two 
initial BCGs cease to dominate their subgroups and one of them comes 
to dominate the merged cluster 



5.1.2. Field subtraction 

To correctly model the field counts, we first need to understand 
the mean surface number density of radio sources in the sky as a 
function of flux. For this we bin the FIRST and NVSS catalogs in 
their entirety in logarithmic flux bins to derive an estimate of the 
number of sources per flux interval and solid angle, d ^ N dsl ■ Due to 
the large sample size of both NVSS and FIRST, this quantity is 
well constrained except at very low and very high flux densities. 
Because it is difficult to find a generally valid fitting function, we 
model d ^ N da using quadratic spline interpolation between bins. 

The number of foreground/background sources, /V« eld , inside 
the search radius of a given cluster is determined by integrating 
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over the field densitjQ as 

where O is the angular area of the searched region, and S min and 
S max are determined from the limiting values of the luminosity 
bin according to Eq. ||3}. 

The total number of cluster sources inside the search radius 

is 

N n = N c ° ant - ld . 

Uncertain ties on N„ are c omputed from Poisson statistics as dis- 
cussed bv lGehrelsl (119861) . 

For reasons of simplicity, the previous discussion has fo- 
cused on source counts in a single cluster field. However, it is 
straightforward to apply the method discussed here to the entire 
stack of clusters by adding all individual N n . 

5.1.3. De-projection into cluster volume 

To convert the cluster source counts N n inside the "line-of-sight 
cylinder" (cf. Fig. |4]l limited by 777-200 to counts N inside 7-200, 
we compute the expected ratio of sources inside the sphere de- 
limited by 7-200 and inside the cylinder. We integrate the spatial 
source density, as a function of de-projected radius, over the two 
regions and take the ratio 

1 

4tt j ¥(H) 3 2 dE 

C ^--T^f--^^ • (6) 

Z=-oop=0 

where *F(S) is the de-projected counterpart of iff(^) (cf. Eq. |4j). 
Note that 3 = rjr^ya, with r now being the de-projected (phys- 
ical) coordinate. The integral over the line-of-sight cylinder (in 
the denominator) is written in cylindrical polar coordinates, with 

Vp 2 + z 2 = 3. 

The de-projected source distribution ^(S) is derived from 
the projected source density distribution i/rtjf), discussed in $3] 
Because no dependence on luminosity could be determined for 
i/f(^), we use the parameters derived from the total sample for 
each combination of cluster and radio source catalogs, as listed 
in Table H 

Because the central Gaussian component of iff(g) in Eq. (HJl 
does not correspond to a physical distribution of sources (cf. 
it is not physically meaningful to de-project it to physical coordi- 
nates. However, because the component is narrow, we can make 
the approximation that all sources belonging to it (i.e. all BCGs) 
have been counted inside the line-of-sight cylinder, provided 77 is 
chosen large enough (77 > 0.1). Thus, no correction in the sense 
of Eq. © is required for the central component. 

We compute a correction factor only for the broader /3- 
model. In order to avoid over-correction where the narrow com- 
ponent is present, we compute from © the relative fraction of 
sources belonging to the broader profile and multiply the correc- 
tion factor by this number to yield a modified correction factor, 
Cn- 



4 We approximate the field as the sum of all regions on the sky cov- 
ered by the radio survey, minus circular regions around the sample clus- 
ters, defined by the respective radii (r 2 oo) estimated from scaling rela- 
tions. 



For the isothermal /3-model, th e physical (de -projected) den- 
sity of the sources is given bv (e.g. lSarazinlll988l) 

Y(S) = fJl + =JJ 2 , (7) 

which is the function we use when computing the integrals in 
Eq. ©. The parameters (/3, £) are unique for each galaxy cluster 
sample, as discussed in |3] Because we are only interested in 
ratios, we set = 1 . 

The cluster volume is defined as a sphere with radius 7-200- 
Thus, the RLF in luminosity bin AL can be expressed as 

(/>(L)AL = Cjy— (8) 

Vsph 

where V sp h is the sum of physical cluster volumes (within 7-200X 
and N n is the total number of sources (in all cluster fields) found 
within a radius 7/ 7-200 and having luminosities in the luminosity 
bin AL. 

The uncertainty in the RLF is dominated by Poissonian noise 
in the source counts. We compute the uncertainty by scaling 
from the uncertainty 6N n in N n as 

6((f>(L)AL) _ 6N n 
(/>(L)AL ~ N,,' 

By construction, the RLF is insensitive to the chosen value 
of the projected fractional radius 77 after the correction factor 
Cn has been applied. However, uncertainties in the RLF can be 
minimized by carefully tuning this parameter. Values too close 
to 1 will increase the background counts, and thereby the error, 
since there are few cluster sources close to 7-200- On the other 
hand, too low values of 77 will increase uncertainties due to not 
providing enough cluster sources for good statistics. We find that 
77 = 0.5 provides a good balance between cluster and background 
sources, and use this value in computing the RLF. 

5.2. Confusion 

Confusion becomes important when the typical separation of 
sources is comparable to the angular resolution of the radio sur- 
vey from which the RLF is derived. In particular, a larger beam 
tends to overestimate the RLF at the high luminosity end as sev- 
eral confused low luminosity sources appear as fewer sources 
with higher luminosity. This effect is redshift dependent through 
the conversion of angular to physical distance, which compli- 
cates quantifying a possible redshift evolution in the RLF. 

Comparing the narrow radial distribution of sources (section 
with the resolution of the FIRST survey, it is apparent that 
confusion affects the counts also when using FIRST data. This 
will affect the shape and normalization of the RLF in a way that 
we have no direct way of quantifying. However, our main inter- 
est is to constrain the redshift evolution of the RLF, and for this 
purpose the absolute normalization of the RLF is of less impor- 
tance. Therefore, rather than attempting to correct for confusion, 
we can correct for the relative difference in confusion at different 
redshifts. 

5.2.1 . Degrading the resolution of the radio data 

To address the confusion problem, we thus degrade the resolu- 
tion of the radio source data in an adaptive way to make sure 
that confusion effects are the same at all redshifts (disregarding 
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second order effects such as a number density evolution in the 
central parts of clusters influencing the confusion problem). In 
the following, we discuss how the resolution of the radio source 
data can be degraded. To demonstrate the method, we apply it to 
the FIRST catalog to enable a direct comparison of the RLF de- 
rived from FIRST to that derived from NVSS, before discussing 
how the method is used in removing first-order systematic effects 
of confusion in the redshift evolution of the RLF. 

The full-width half maxima (FWHM) of the FIRST and 
NVSS synthesized beams are 5" and 45", respectively. To what 
extent individual sources can be distinguished in each of the sur- 
veys, however, largely depends on their absolute and relative 
brightnesses. Here we make the simplification that given a reso- 
lution in terms of a beam FWHM, no two sources are resolved 
if their mutual separation (irrespective of flux) is less than the 
FWHM. 

To degrade the resolution of the FIRST or NVSS data, we 
define a beam with a FWHM larger than that of the original data. 
In a given field, the position and flux density of each source are 
recorded. The brightest source is located, and all sources within 
a radius equal to the FWHM around this source are combined 
into one source. The position of the new source is taken as a 
flux-weighted average of its parts, and the total flux density is 
the sum of the integrated flux densities of the parts. The next 
brightest source (excluding all sources that have already been 
considered) is then located, and the procedure is repeated until 
no two sources are separated by an angular distance less than the 
FWHM. Note that this method is limited by the use of source 
catalogs in the sense that sources fainter than a chosen limit will 
not be considered. Here we use the completeness limit of the 
two source catalogs to determine which sources to consider for 
summing. 

As an improved method, one could instead use the raw maps 
of cluster fields from the FIRST survey and degrade them to 
the NVSS resolution before extracting sources. However, this 
requires a robust source extraction in both maps and is quite 
complicated for faint sources near the completeness limit of the 
surveys. Moreover, neither such a method nor the method de- 
scribed above can resolve the inherent problem that some ex- 
tended emission recovered by NVSS is resolved out by FIRST. 

5.2.2. NVSS and FIRST source counts 

A simple way to visualize the effect of confusion on the RLF is 
to compare the luminosity function as constructed from NVSS 
with that constructed from FIRST, using the same complete- 
ness cutoff in both surveys. For this comparison we select clus- 
ter fields in a narrow redshift range, 0.1 < z < 0.16, from the 
maxBCG catalog to ensure a sufficiently large statistical sample. 
This sub-sample of the maxBCG catalog is used only for the 
purpose of directly comparing luminosity functions constructed 
from FIRST and NVSS in this section. It contains 2341 clusters 
with M 20 o > 5 x 1O 13 M which are covered by NVSS and FIRST. 
As indicated in Fig. [5] computing the RLF directly from the two 
radio source samples yields inconsistent results. This is expected 
due to confusion effects. 

We attempt to re-create the NVSS based luminosity function 
from the FIRST data by degrading the resolution as described 
above. Note that we do not expect a perfect agreement since the 
sensitivity to spatial frequencies of FIRST is very different from 
that of NVSS. 

Because both surveys have many sources separated by dis- 
tances much smaller than the respective resolutions (due to the 
specific source extraction methods applied to the synthesized im- 



ages), we impose a strict lower limit of 45" on the separation 
of sources in both surveys for this comparison. The results are 
shown in Fig. [5] Even using this simple method, the RLF from 




23.0 23.5 24.0 24.5 25.0 25.5 26.0 26.5 

Log[P (W Hz" 1 )] 



Fig. 5. Change in the RLF resulting from degrading the reso- 
lution of both the FIRST and NVSS catalogs to a minimum 
source separation of 45" (the FWHM of the synthesized beam 
of NVSS). The shaded regions indicate the RLF computed from 
the non-degraded NVSS (grey) and FIRST (hashed) data. The er- 
ror bars indicate the degraded versions; here the counts derived 
from the FIRST catalog (blue, thin error bars) are in approxi- 
mate agreement with the NVSS counts (green, thick error bars), 
despite the large difference in synthesized beam FWHM of the 
original data. 



the degraded FIRST sample is in good agreement with that de- 
rived from the NVSS counts. The RLF constructed from the de- 
graded FIRST data is systematically lower than that constructed 
from the degraded NVSS data. This is as expected considering 
the different sensitivities to spatial frequencies of the two sur- 
veys - extended emission recovered by NVSS is resolved out by 
FIRST, causing the amplitude of the RLF to drop. Note, how- 
ever, that the difference in the luminosity function at different 
redshifts caused by this effect will in fact be much smaller than 
the residual difference seen in Fig. [5] since our relatively small 
redshift range in the maxBCG sample results in a much smaller 
relative difference in physical scales, as compared to the relative 
difference between the FIRST and NVSS synthesized beams. 

5.2.3. Accounting for confusion 

Given the above model, it is possible to adaptively introduce 
confusion into our radio source sample to minimize systematic 
effects in determining the redshift dependence of the RLF. In a 
given cluster sample, the clusters with the greatest angular di- 
ameter distances will be the most affected by confusion. Thus, 
we define a nominal "confusion distance", in physical units, by 
converting the radio source survey resolution (45" for NVSS or 
5" for FIRST) to a physical distance at the redshift of the ob- 
ject with the greatest angular diameter distance in the sample. 
Then, for every object in the sample, this physical distance is 
converted back to angular units using the redshift of the object, 
and the resulting angular scale is used to degrade the resolution 
of this particular cluster field using the method described above. 
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To get an idea of the typical level of confusion, consider the 
maximum redshift, z = 0.3, of the maxBCG catalog. The FIRST 
resolution of 5" corresponds to a physical distance of 22 kpc at 
this redshift. We can translate this physical distance, the "confu- 
sion distance", into an angular distance at any lower redshift; for 
example, at the low redshift limit of z = 0.1 we obtain an angular 
scale of 12". Using this as a limit of resolution, the FIRST source 
counts are reduced by approximately 25% in fields defined by 
7"200> while at z = 0.3 the FIRST source counts are reduced by 
around 10%Q 

Note that we also have to adapt the background counts, A^ eld , 
when degrading the resolution. Because it is time consuming 
to re-compute d $ N da for each angular resolution used, we com- 
pute the field counts in several degraded versions (as described 
in 35 .2. It of the entire NVSS and FIRST catalogs, using incre- 
ments of 5", and interpolate between these results to derive esti- 
mates of A^ eld individually for all cluster fields. 

5.3. Modeling the luminosity function 
5.3.1. Parametric model 

To allow for a quantitative estimate of the redshift evolution of 
the RLF, we follow LM07 and use the parameterization 



log 4> = y - 



2 ( log.L-* 

D + 



2\ 



1/2 



1.5 log L. 



(9) 



This fitting function was used by ICondon et al.l (l2002h to fit 
the field RLF using the combined contributions from radio-loud 
AGN and star-forming (SF) galaxies. LM07 used the two sets 
of values of (b, x, w) found by Condon et al. to make possible a 
separation of the normalizations (y) of the radio-loud AGN and 
SF components in their cluster RLF at low redshift. 

It is beyond the scope of this paper to discuss the physical in- 
terpretation of Eq. (|9} as our main concern is the r edshift evolu- 
tion of the luminosity and number density. Refer to lCondon et al.l 
(2002) and references therein for details on this hyperbolic fit- 
ting form of the RLF. 

In contrast to the sample studied by LM07, our main samples 
do not include systems with z < 0.1. This selection, combined 
with our flux limit of 5 mJy, effectively ensures that we are only 
sensitive to the high-luminosity radio-loud AGN population of 
the RLF. Thus, a one-component model as described by Eq. <j9j 
is sufficient for the present purpose. 

5.3.2. Redshift evolution 



iMachalski & G odlowskil (l2000l) quantified the redshift evolution 
of the RLF under the assumpti on that th e overa ll shape remains 
constant, as first suggested by Condon] (Il984l) . Under this as- 
sumption, there can only be changes in overall luminosity and 
overall number density. The redshift dependence can be written 
as 

<f>(L,z) = g(z)c/>[Lf(z),z*0], (10) 

where g(z) quantifies the number density evolution, correspond- 
ing to a vertical shift AY in tp(L), and f(z) represents luminosity 
evolution, corresponding to a horizontal shift AX: 

log [0(log L, z)] = log [0(log L + AX, z ~ 0)] + AY. (11) 



5 The counts are reduced also at the resolution limit because, contrary 
to the FIRST catalog, we enforce a strict limit of 5" as the smallest 
angular distance between two distinct sources. 



The vertical and horizontal shifts can be fitted to yield f(z) and 
g(z). Shape preservation of the RLF in (log L, log <f>) space im- 
plies that number density and luminosity scale as powers of red- 
shift. For the luminosity, 



L = Lq 



1 +z 

1 +Z0 



(12) 



where ai is the power law index, and correspondingly for the 
number density 

(1+7 
T+I0/ ■ (13) 

with the power law index a^. 

To constrain the redshift evolution in the RLF, we bin our 
samples by redshift, and approximate the redshift in each bin by 
the median of all cluster redshifts in the bin. We then apply 
shifted in <p and L according to (fTTT i. to constrain the power law 
indices and a^. Simultaneously, the four parameters of (|9) are 
constrained. 

We let the parameter y in (0 represent the normalization of 
the RLF at z = 0. We bin the data by redshift and construct the 
RLF separately for each redshift bin. The thus constructed lu- 
minosity function data are then used simultaneously to constrain 
y, b, x, w, a,j, and The fitting is carried out by minimization of 
the^ 2 statistic, 

k lv \ 2 



X 



S 

i=l 



(14) 



where X-, is the data point with index i (the number k of data 
points is the total number of luminosity bins in all redshift bins), 
Hi is the corresponding model value given a set of parameters 
(y, b, x, w, a,/, and c?i), and <x ; is the uncertainty on data point ; 
estimated from Poisson statistics. 

Uncertainties on the fitted parameters are estimated using a 
simple method. Given our data points, binned in redshift and 
in luminosity, we vary the data according to Poisson statistics 
within their estimated uncertainties, and carry out the x 2 fit 
again. This process is repeated 10.000 times. As expected, the 
mean of the distribution of a fitted parameter in any of the fits 
we carry out ( {35.4. 21 ) corresponds to its best-fit value. The stan- 
dard deviation of the distribution of a fitted parameter is used as 
an estimate of the lcr uncertainty of this parameter. 

5.3.3. Mass dependence 

A possible mass dependence in the RLF can be investigated us- 
ing the same method as described above. By scaling the num- 
ber density and luminosity we can model the dependence to a 
first approximation as power laws of the form M y my Analogously 
with the discussion in 35.3.21 we make the assumption of shape 
preservation of the RLF and model the luminosity and number 
density dependence on mass as (cf. Eqs. (fT2l) and ( TP3| )) 



L ~ (M 20 o) ri ; 

4>~(M 2m y*. 



(15) 



We use both the X-ray sample and the optical sample to con- 
strain the mass dependence. Because two new parameters are 
introduced into the model, it is not possible to simultaneously 
constrain the mass dependence and the redshift evolution, par- 
ticularly given the limited size of the X-ray sample. For this rea- 
son we use a sub-sample of the X-ray selected clusters in a low 
redshift range of 0.05 < z < 0.12 (see Table[2]i. In this range, we 
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Fig. 6. RLF determined from the X-ray subsample in the red- 
shift range 0.05 < z < 0.12, divided into two mass bins. Arrows 
indicate lcr upper limits from Poisson statistics. 
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Fig. 7. RLF determined from the optical sample, divided into 
three mass bins. 



divide the sample into two mass bins, where care is taken that the 
lower mass limit of the lower mass bin is above the completeness 
limit of the X-ray sample at z — 0.12. For the optical sample, we 
divide the complete sample into three mass bins. In the next sec- 
tion it is shown that in this range, our data are consistent with no 
mass dependence of the RLF. 



Table 7. Mass dependence of the RLF, parameterized in terms of 
power laws of M200 in luminosity and number density (Eq. (15[). 



Cluster 


Source 


Prior 


y* 


JL 


sample 


catalog 








X-ray 


FIRST 




0.25+1.01 


-0.46±0.55 


maxBCG 


FIRST 




-0.36+0.31 


0.38±0.65 


X-ray 


FIRST 


7> = 


(0.0) 


0.13±0.36 


maxBCG 


FIRST 


7* = 


(0.0) 


0.27±0.28 


X-ray 


FIRST 


7l = 


0.090±0.16 


(0.0) 


maxBCG 


FIRST 


7l = 


0.046±0.15 


(0.0) 
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Fig. 8. Redshift evolution in the RLF determined from the op- 
tical (maxBCG) sample, without priors. The data are shown as 
error bars and the best simultaneous fit to luminosity and num- 
ber count evolution is represented by lines corresponding to 
the redshift bins: 0.1 < z < 0.17 (diamonds and solid line), 
0.17 < z < 0.24 (triangles and dashed line), and 0.24 < z < 0.3 
(squares and dotted line). 



Thus, contrary to LM07, who find that the amplitude of 
the RLF for low-mass clusters is slightly larger at the low- 
luminosity end, we find no evidence of a mass dependence in 
the RLF given our flux and redshift cut-offs. As indicated in 
Figure [6] the results are consistent with the most luminous ra- 
dio sources only being found in high-mass clusters, as was also 
found by LM07. However, due to the relatively small sample 
size, our data provides no conclusive evidence of this. 

As there is no statistically significant indication of a mass 
dependence from either cluster sample, we proceed to fit for a 
pure redshift evolution in the next subsection. 



5.4. Results 

5.4.1. Mass dependence 

To separate a possible mass dependence in the RLF from a pure 
redshift dependence, we construct RLF binned by mass as shown 
in Figures [6] and Q We then use a least-squares statistic to find 
the best-fit vertical and horizontal shifts of the RLF, fixing the 
values of b, w, x as described in 35.3.2I to interpret the shifts 
in terms of a luminosity- and/or number density dependence on 
cluster mass. The results for both our cluster catalogs are given 
in Table |7l using different priors on the power laws. Note that all 
fits are consistent with — and ji — at an approximate 1 <x 
level. 



5.4.2. Redshift evolution 

Figure [8] shows the RLF computed from the maxBCG sample, 
using FIRST radio sources, in three redshift bins chosen such 
that there are approximately the same numbers of clusters in 
each bin (0.1 < z < 0.17,0.17 < z < 0. 24 and 0.24 < z < 0.3). 
We apply the method outlined in jH5.3.2l to constrain the redshift 
evolution of the RLF by a simultaneous fit to the RLF amplitude 
(parameterized by y at z = 0), the RLF shape (parameterized by 
b, x and w), and the power law evolution (parameterized by a^, 
and ai according to Eqs. (fT2l and (fTTt). 

We perform a number of fits to the data, with priors as listed 
in Table [8] Although a simultaneous fit to both and yields 
the best fit (in the sense of the reduced x 1 parameter, x^,, be- 
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Table 8. Redshift evolution in the RLF, quantified by simultaneous fits to the RLF amplitude, shape and evolution. The goodness- 
of-fit is indicated by the reduced chi-squared parameter, xLj- See section l531 for a description of the parameters. 



Cluster 


Source 


Priors 


y 


b 


X 


w 


a> 


OIL 


'i 


sample 


catalog 
















maxBCG 


FIRST 




36.38+1.02 


1.05+0.73 


24.53+0.18 


0.66+0.13 


-2.46+1.58 


6.20+1.76 


1.07 


maxBCG 


FIRST 


a$ = 


36.34±0.92 


0.91+0.81 


24.87+0.14 


0.72+0.21 


(0.0) 


3.99+1.24 


1.19 


maxBCG 


FIRST 


a L = 


36.74±0.89 


1.01+0.55 


25.11+0.11 


0.71+0.19 


1.03+1.14 


(0.0) 


2.25 


X-ray 


FIRST 


(a) 


36.19+0.19 


(1.05) 


(24.53) 


(0.66) 


0.76+1.86 


8.12+2.67 


0.94 


X-ray 


FIRST 


(a); a = 


36.26+0.10 


(1.05) 


(24.53) 


(0.66) 


(0.0) 


8.19+2.66 


0.89 


X-ray 


FIRST 


(a); a L = 


35.89+0.18 


(1.05) 


(24.53) 


(0.66) 


9.40+1.85 


(0.0) 


10.48 



(a) The shape parameters b, x and w were fixed to the results of the maxBCG/FIRST analysis with no priors. 



ing the closest to unity), the data are consistent with no number 
density evolution. Note that the best-fit value a$ = -1.38 im- 
plies fewer radio-loud AGN at higher redshift within the range 
0.1 < z < 0.3. Though not unphysical, such a negative evolution 
is unlikely, considering the evidence of increase d AGN activity 
in the field population (Dunl op & Peacocklll990h . and the en- 
hanced AGN fraction within clusters as seen from X -ray obser- 
vations (East man et al.l 120071: iGala metz et aTl l2009t) Therefore, 
we carry out an additional fit in which we fix the number density 
power law to zero and fit for a pure luminosity evolution. The 
result is a positive luminosity evolution with ai = 3.99 + 1.24. 
While consistent with a pure luminosity evolution, the data are 
inconsistent with a pure number density evolution, as indicated 
by the^ ed of the third column in Table [8] 

The last three rows of Table [8] indicate the results of fitting 
the evolution parameters to the X-ray derived RLF, also binned 
by redshift (redshift bins: 0.1 < z < 0.2, 0.2 < z < 0.45 and 
z > 0.45). Although the X-ray sample extends to greater red- 
shifts than the maxBCG sample, it has less leverage on the red- 
shift evolution because of its limited size. For this reason, we fit 
only for y, a$ and a^. We keep the shape parameters fixed to 
the values found in the maxBCG/FIRST fit, as this is the fit with 
the smallest uncertainties. Because the smallest xLa is measured 
for the case of no priors (as expected) in the maxBCG/FIRST 
case, we use the results of this fit for fixing the shape parame- 
ters. Again, we also carry out fits to pure luminosity evolution 
and pure number density evolution. At one standard deviation, 
the results from the X-ray sample are consistent with those of 
the optical maxBCG sample. Again, a pure number density evo- 
lution provides a poor fit to the data. 

We note that in all cases the fit is especially sensitive to the 
parameter b, which is explained by its quadratic dependence in 
Eq. @. 

5.5. Systematic uncertainties 

To estimate systematic uncertainties in the luminosity and num- 
ber density evolution of the RLF, we discuss below four separate 
sources of error. 

The estimated systematic uncertainties for the maxBCG and 
X-ray samples are summarized in Table |9] For each uncertainty 
estimate, we keep all parameters but the RLF normalization and 
the evolution parameter in question (a^ or ai) fixed, while car- 
rying out the complete RLF analysis with posterior assumptions 
altered as described above. To estimate total systematic effects 
on and a/,, we add all relevant components in quadrature. 
We note that for this study, the statistical errors on the fitted pa- 
rameters (see Table|8) still dominate over the systematics. 



Table 9. Summary of systematic effects on the number density 
and luminosity evolution of the RLF. 



Systematic 


effect on a$ 


effect on a L 


L bo i - T scaling 


±0.30 




L x - n™° scaling^ 
fc-correction 


±0.80 

+0.15 
-0 51 


+0.17 
-0 16 


Totals" 


maxBCG 


+0.87 


+0.17 


X-ray 


-M 

-n.59 


M 

-n if, 



" only for maxBCG clusters 
b added in quadrature 



5.5.1. Scaling relations 

Uncertainties in the z-dependence of the scaling relations used 
to derive cluster masses affects the derived number density evo- 
lution through the indirect determination of cluster volumes. For 
the X-ray sample, we estimate this systematic by replacing the 
self-similar evolution of the Lboi ~ T rela t ion by (i) the strong 
evolution measured bv lKotov & Vikhlininl (l2005) and (ii) by as- 
suming no evolution of the Lboi - T relation. Re-computing the 
M200 for both cases results in a change in the number density 
evolution parameter a$ of +0.30 when re-computing the RLF 
for the X-ray sample. For the optical sample, we additionally 
need to take the uncertainty in the Lx - n^ 00 relation (Eq. (JTJ) 
into account. We estimate the systematic effects by using the ex- 
treme values (at lcr c onfidence) of the redshift dependence from 
iRvkoff et a l. (2008b) to re-compute M200. resulting in an effect 
on of ±0.80. 

5.5.2. ^-correction 

Uncertainties in the ^-correction used in converting from radio 
source flux to luminosity will mainly affect the luminosity evolu- 
tion. It can also shift the number density evolution by assigning 
the wrong luminosity bins to sources. 

As extreme values of the ensemble average spectral slope, 
we use the 25th and 75th percentiles (a - 0.51 and a — .92, re- 
spectively) of the slope distribution of lCoble et al.l(l2007l) . We re- 
compute the maxBCG RLF using these extreme spectral slopes 
for all sources. 

5.5.3. Source counts 

Incorrect radio source counts can be caused either by counting 
extended or complex sources (e.g. AGN with radio lobes) as sev- 
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eral separate sources, or by not resolving unrelated sources. This 
affects the number density and luminosity evolution only if the 
effect is varying with redshift; otherwise a mere offset in the nor- 
malization of the RLF is introduced. 

A complete treatment of these effects would require the care- 
ful visual inspection of the radio maps in each cluster field, 
which is not feasible for a large sample such as maxBCG. 
Instead we carried out a visual inspection on parts of the sample 
to estimate the possible effects. We visually inspected FIRST 
maps in the 500 fields corresponding to the clusters with the 
lowest redshifts, and in the 500 fields corresponding to the clus- 
ters with the highest redshifts of the maxBCG sample, taking 
weighted sums of the flux densities of FIRST catalog entries 
deemed to be different components or lobes of the same sources. 
Ambiguous cases (less than 2% of the fields studied) were re- 
moved. Comparing the RLFs constructed from the modified 
FIRST data in the two sub-samples, we noted a drop in normal- 
ization of the RLF by about 5% in both the high-z and the low-z 
subsamples. There was no indication of a significant change to 
the evolution of the RLF within our statistical errors. Thus, we 
conclude that our adaptive accounting of confusion with redshift 
(section l5T2l is sufficient to deal with this systematic, at least for 
the purpose of constraining the redshift evolution of the RLF. 

5.5.4. Radio source flux 

Systematic offsets in radio source flux densities are expected to 
be a problem mainly for the FIRST survey, where significant 
amounts of flux of extended sources can be resolved out. Again, 
this problem only affects the normalization of the RLF, unless 
there is a redshift dependence in the amount of resolved-out flux. 
Since the results from the NVSS and FIRST radio data for the 
redshift evolution of the RLF from the maxBCG sample are con- 
sistent within statistical errors, we do not pursue this point fur- 
ther. 



5.6. Comparison of the RLF with previous findings 

To compare our findings with previously published results, we 
consider the maxBCG/FIRST data in our lowest redshift bin, 
which compares well to the low -redshift samples of clusters used 
in previous works. Figure |9] shows a comparis o n of our low- 
redshi ft R LF to the results oflLin & Mohrl (120071) . iReddv & Yunl 
(12004 and lMassardi & De Zottil (120041) . 

We note that the normalization of the RLF is consistent with 
the results of LM07, although the RLF from our optically se- 
lected sample does not extend as high in luminosity (Log(L) > 
26 WHz 1 ). The latter is likely a result of the fact that the op- 
tically selected sample contains many more low-mass clusters, 
and that the more luminous radio sources tend to reside in high- 
mass clusters (cf. {J5AJ}. 

The RLF of Massard i & De Zottil d2004t) has a significantly 
lower normalization than what we find, by approximately a fac- 
tor of two. This can be expected due to the different definition of 
cluster volumes (using a constant cluster radius across the com- 
plete sample) used in the former work. 

iBranchesi et alJ (2006) used a definition of the RLF as the 
number of radio sources per cluster rather than averaged over 
a volume. A dire ct comparison w i th that work would require a 
re-analysis of the IBranchesi et al.l (120061) sample, which would 
be only of moderate interest considering that the normalization 
would rely on different assumptions on cluster volume and the 
sample is selected differently from our maxBCG sample. Our 
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Fig. 9. Comparison of our RLF derived from the 
maxBCG/FIRST data in the lowest redshift bin (0.1 < z < 0.17, 
Zmedian = 0.13; black diamonds) with other rele vant volume- 
averag ed radio luminosity functions. The RLF of Lin & Mohr 
(2007) is indicated by shaded regions, while the c orre- 
sponding luminosity funct ions of IReddv & Yunl d2004l) and 
iMassardi & De Zottil d2004l) are shown as green triangles and 
red squares, respectively. 



main objective here is to show the approximate agreement be- 
tween different local luminosity functions and to stress the dif- 
ficulty in obtaining a robust normalization of the RLF; although 
the normalization depends strongly on the choice of cluster vol- 
ume and is prone to systematic effects due to a poor understand- 
ing of this volume, this does not strongly affect our main conclu- 
sion, which is concerned with the redshift evolution of the RLF 
independently of its overall normalization. 

Early studies of the luminosity function of optically selected 
QSOs at z < 2.2 suggested a pure lum inosity evolution with L 



(1 + z )(3-s±0.3) in the field p p U i at i on dBovle et alJl987lll988h. A 



later s tudy of extragalactic radio sources by Dunlop & Peacock 
(1990) at 2.7 GHz came to similar conclusions. 

The pure luminosity evolution of the volume averaged RLF 
found in this study is approximately consistent with more recent 
findings in the field population of radio galaxies at low and in- 
terme diate redshift (Machals ki & Godl owski 2000; B rown et"aTI 
l200lh . although our best-fit evolutionary model with L ~ 
(1 + z )(6- 2Q±i-76) i s steeper t han both the L ~ (1 + z) (4±1) 
found bv iBrown et al.l d200ll) and the L ~ (1 + z) (3±1) which 
Machalski & Godlowski (2000) found to be consistent with their 
data. 



6. Summary and conclusions 

We used the maxBCG optical sample of galaxy clusters and 
a composite X-ray sample to construct the volume averaged 
radio luminosity function (RLF) in galaxy clusters by cross- 
correlating cluster positions with radio point source positions 
from the FIRST and NVSS survey catalogs. Background and 
foreground counts were corrected for, and variable confusion 
with redshift was accounted for. We investigated the radial 
source density distribution of radio sources associated with clus- 
ters, and correlated the luminosity of the brightest radio source 
with cluster mass. 

To fit the radial number density distributions, we used 
the functional form of the isothermal y6-model. Combining the 



14 



M. W. Sommer et al.: Galaxy cluster radio luminosity function 



maxBCG cluster sample with the FIRST catalog, we found an 
additional narrow component which we identify with radio de- 
tections of the brightest cluster galaxies (BCGs). All combina- 
tions of cluster sample (maxBCG or X-ray) and radio source 
sample (NVSS or FIRST) yield similar values of the power law 
index /3. The core radius in the maxBCG/FIRST radial distri- 
bution of sources compares well to the distribution of sources 
in the X-ray sample. As it was not possible to identify the nar- 
row BCG component from the NVSS survey, we found a much 
smaller core radius in the maxBCG/NVSS radial distribution. 

Our derived radial distributions are narrower than what 
has been found in earlie r studies (e.g iReddv & Yunl [2004; 
Mass ardi & De Zottill2004l) . even when disregarding the central 
component associated with brightest cluster galaxies (BCGs). 
A plausible explanation is that we are constructing the RLF at 
higher redshifts than previous studies and thus, given our flux 
limit of 5 mJy, we are sensitive only to radio-loud AGN. At lower 
redshift there is a mixing with star forming galaxies, which have 
a less centralized distribution. 

Unlike LM07, we do not find any evidence that bright radio 
sources have a radial source density distribution different from 
that of faint sources. Again, a likely explanation is that we are 
studying different populations of radio sources through the red- 
shift selection. 

We found that the luminosity of the most radio-luminous 
source within 50 kpc from the cluster center scales with clus- 
ter mass following a power law with slope 0.31+0.12 in 
the maxBCG samp le. This is consistent with the results of 
iLin & Mohri £2004) as well as with the results from our X-ray 
sample, although the latter is also consistent with no correlation. 

We find the RLFs constructed from the optical and X-ray 
samples of galaxy clusters to be in approximate agreement. The 
RLF from the optical maxBCG sample is systematically lower at 
luminosities L > 3 x 10 25 W Hz 1 . This is likely a result of many 
more low-mass systems being present in the optical sample. 

We provide the first evidence for a luminosity evolu- 
tion of the volume-averaged RLF in clusters of galaxies. The 
maxBCG/FIRST data are consistent with a pure luminosity evo- 
lution, with power scaling with redshift as L ~ (1 + zf L , where 

= 6.20*Hgj|? (statistical followed by systematic uncertain- 
ties). There is no indication of a mass dependence in the RLF 
from the present data. However, the results from the X-ray sam- 
ple are consistent with the findings of LM07, that the most lumi- 
nous radio sources reside in massive clusters. This is further cor- 
roborated by the fact that the RLF constructed from the maxBCG 
sample (which contains a smaller fraction of high-mass systems 
than both our X-ray sample and the sample of LM07) is steeper 
at higher luminosities. 

Below P ~ 10 25 WHz~', the data are consistent with no 
mass dependence in the RLF, as shown by constructing the RLF 
in the low-redshift X-ray sample and binning by mass. Although 
we have found that massive clusters have more luminous BCGs, 
this effect is counteracted in the RLF by the fact that these mas- 
sive systems also have more volume. Both the X-ray sample and 
the maxBCG sample yield results consistent with a pure lumi- 
nosity evolution of the RLF. In addition, the derived power laws 
are comparable, although the samples cover very different red- 
shift ranges. 
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